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ABSTRACT 

We measure the zero point of the Stanek (1996) extinction map by comparing 
the observed (V — K ) colors of 20 RR Lyrae stars (type ab) found in the 
MACHO survey with their intrinsic (V — K ) 0 colors as a function of period 
as determined from nearby RR Lyrae stars. We find that the zero point of 
the Stanek map should be changed by A Ay = —0.11 ± 0.05 mag, in excellent 
agreement with the recent measurement of Gould, Popowski & Terndrup (1997) 
using K giants. 

Subject Headings: dust, extinction - Galaxy: general - stars: variables: 
other 
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1. Introduction 

Baade’s Window, (£,b) ~ (1°, —4°), has been an important laboratory for the study 
of bulge populations because of its relatively low extinction {A v ~ 1.5) and the presence 
of the globular cluster NGC 6522 which provided an early opportunity to measure that 
extinction. Baade’s Window has been a major focus of microlensing surveys by OGLE 
(Udalski et al. 1993), MACHO (Alcock et al. 1995), and EROS (Ansari et al. 1996). Among 
the many non-microlensing applications of these surveys is Stanek’s (1996) construction 
of a detailed extinction map with 30" resolution over the 40' square OGLE field. Stanek 
(1996) applied the method of Wozniak & Stanek (1996), in essence measuring the mean 
apparent magnitude, (V), and color, (V — /), of clump giants as a function of position. He 
then inferred from these the differential total and selective extinctions, Ay and E(V — /). 
The estimated relative error in Ay of this map is 0.1 mag for individual resolution elements, 
but there is an overall uncertainty of 0.2 mag in the zero point. 

For many applications, such as the interpretation of color-magnitude diagrams of bulge 
held stars and of the cluster NGC 6522 or the measurement of distances using RR Lyraes 
or other tracers, the determination of the zero point is crucial. 


2. Previous Work 

2.1. E(B — V) 

Historically, the extinction toward Baade’s Window has been found by measuring 
the selective extinction E{B — V) , and then multiplying by an assumed ratio of total to 
selective extinction Ry = Ay/E{B — V ) (Arp 1965; van den Bergh 1971; Walker & Mack 
1986; Terndrup & Walker 1994). There are several major disadvantages to this approach. 
First, Ry ~ 3 is rather large, and the statistical error in E{B — V) (usually estimated to 
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be > 0.03) is multiplied by this factor when estimating the error in Ay. Second, Ry varies 
along different lines of sight, so for any particular line of sight for which it is not actually 
measured (e.g., Baade’s Window) the precision of the estimate is no better than 7%. Hence, 
the statistical error alone for Ay is more than 0.12 mag. Finally, there are systematic errors 
arising from uncertainties in the intrinsic ( B — V) colors of stars used to estimate E(B — V). 
While the intrinsic color of extremely hot stars (in the Raleigh-Jeans limit) is known from 
fundamental physics, there are no known stars at such extreme temperatures lying beyond 
the dust column in this direction. Hence, one must use cooler stars whose (B — V ) colors 
are sensitive functions of temperature, metallicity, and perhaps other factors. The standard 
approach is to find local analogs of the program stars and directly measure their colors, but 
systematic errors may arise from any unrecognized differences between these two groups of 
stars. As always, it is difficult to determine the size of the systematic errors, but one can 
gain a sense of their magnitude by comparing the E(B — V)q = 0.45 ± 0.04 derived by van 
den Bergh (1971) from three different methods based on cool stars (K and M giants) with 
the E(B — V) 0 = 0.60 ± 0.03 derived by Walker & Mack (1986) using relatively hot stars 
(RRab Lyraes). Here the subscript “0” means “reduced to zero color” using the prescription 
of Dean, Warren, & Cousins (1978). 


2 . 2 . E(V - K) 

Terndrup, Sadler, & Rich (1995) pioneered a radically different approach. They 
measured E(V — K) = 1.23 ± 0.08 for a sub-region of Baade’s Window (“Blanco 
region A”) with relatively uniform extinction by comparing the H/3 A4861 (Faber et 
al. 1985) index as a function of (V — K) color to that observed for bright K giants 
in the solar neighborhood. They then inferred (but did not explicitly write down), 
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A v = E(V — K)/( 1 — a) = 1.38 ± 0.09, where 

a 

a = —— = 0.112 ± 0.002, (1) 

Ay 

is the ratio of K to V extinction (Rieke & Lebofsky 1985). While this approach is formally 
identical to the previous one (measurement of a selective extinction and conversion to a total 
extinction), it is potentially more accurate than using E(B — V) because the extrapolation 
to total visual extinction is small (a factor 1.12 vs. 3), and therefore the error in Ay is only 
slightly bigger than the error in E(V — K). 

Stanek (1996) set the zero point of his extinction map by forcing it to reproduce the 
results of Terndrup et al. (1995) over their sub-region. However, as discussed in some 
detail by Gould, Popowski, & Terndrup (1997), Stanek (1996) incorporated inconsistent 
assumptions in setting the zero point. Gould et al. (1997) therefore estimated a “naive” 
correction to the Stanek (1997) map of A A v = —0.09 ± 0.09. 


However, Gould et al. (1997) also recognized that the very existence of the Stanek 
(1996) differential map makes possible a much more accurate determination of the zero 
point. One can now use the individual differential extinctions Hy? nek from the Stanek (1996) 
map for an ensemble of stars i = 1 ...n to make individual estimates for the correction to the 
Stanek (1996) map, 


A A Vji — 


C V - K)i - ('V - K) 0 ,i 


a 


a Stanek 

^ V,i 


( 2 ) 


Here (V — K)i is the observed color of the star and (V — K) 0) j is its predicted dereddened 
color. Gould et al. (1997) applied this method to a sample of n = 206 K giants from 
Terndrup et al. (1995) and derived 


AAy = (A Ay ti ) = -0.10 ± 0.06 (K Giants). (3) 

Here we apply equation (|2j) to a sample of n = 20 RR Lyrae stars (type ab) found by the 
MACHO collaboration. RR Lyrae stars are substantially hotter than K giants and the 
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method by which we infer their dereddened (V — K ) 0 colors is substantially different than 
that used by Gould et ah (1997). Therefore, this new determination yields an important 
check on any unrecognized systematic effects which may have affected previous results. 

3. Calibration 

Figure 1 shows the dereddened (V — K) 0 colors for 20 RR Lyraes plotted against the 
fundamentalized period, P 0 , taken from Table 9 of Jones et al. (1992). For the 17 RRab’s 
(circles), this quantity is the same as the observed period, while for the 3 RRc’s ( crosses ) it 
is inferred from the observed first-overtone period P\. The solid line is the best fit to the 
RRab’s and is given by 

(V - K) o = 1.046 ± 0.013 + (1.245 ± 0.144)(logP 0 + 0.29), (4) 

where the expression gives the fit in a form with uncorrelated errors. In more 
traditional format, this becomes (V — K )o = 1.407 + 1.245 log Pq. It is also 
possible to include [Fe/H] as an independent variable, in which case the best fit is 
(V — K) o = 1.650 + 1.7971ogP 0 + 0.084[Fe/H], However, as we show below, there is little 
advantage to this. 

It is clear from Figure 1 that the three RRc’s do not follow the same relation as the 
RRab’s. Moreover, there are not enough RRc’s to determine from the data the relation 
that is appropriate for them. We therefore exclude RRc’s from further consideration. 

4. Data 

The data are drawn from four sources and are presented in Table 1. Columns 1 and 2 
show the positions of each star in 2000 coordinates. Column 3 gives the periods and column 



4 gives the apparent V magnitudes, both from MACHO data (see §4.1). Column 5 gives 
the apparent K magnitudes from Table 1 of Carney et al. (1995) (see §4.2). Column 6 gives 
the metallicity from Walker & Terndrup (1991). Column 7 gives the visual extinction from 
the Stanek (1996) map. Columns 8 and 9 gives the actual (V — K ) (from columns 4,5 and 
7) and the (V — K ) predicted on the basis of equation (|4|) using columns 3 and 7. 


4.1. V Band Photometry 

The apparent V magnitudes listed in Table 1 are each based upon of order 100 
observations taken by MACHO in each of two passbands Bm and Rm, during the 1993 
bulge season. For each star, we find the magnitudes at mean flux, (Bm) and (Rm), by 
taking the average flux (after excluding bad data points) for all observations, and then 
converting to a magnitude. Since the dispersion of a typical RRab light curve is ~ 0.3 mag, 
this procedure produces a random error of ~ 0.03 mag (relative to perfect coverage of the 
light curve) which is small by comparison with other errors. We convert these instrumental 
magnitudes to standard Johnson V band using V = 23.67 + 1.0026F>m — 0.156(Rm — Rm) 
(Alves 1997). We apply this conversion to the magnitudes at mean flux, (B M ) and (Rm)- 

As a check on the zero point of this conversion, we compare in Figure 2 the V 
magnitudes at mean flux derived from MACHO data with those derived by Olech (1997) 
from OGLE data for 50 RRab’s in common between the two data sets. Since Olech (1997) 
reports mean magnitudes, not magnitudes at mean flux, we first correct the OGLE values 
by the difference between the two quantities as determined from the MACHO data. The 
offset of OGLE relative to MACHO is 0.006 ± 0.019 mag, with a scatter of 0.13 mag in the 
difference. That is, the zero points are in excellent agreement. We note that Alcock et al. 
(1997) found a scatter of 0.12 mag in the mean mags determined from MACHO photometry 
of individual bulge RR Lyrae stars by comparing 44 stars detected in overlapping fields. 
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Hence, the scatter in the MACHO/OGLE difference can largely be accounted for by the 
MACHO scatter. 


4.2. K Band Photometry 

The apparent K magnitudes listed in Table 1 are taken from Carney et al. (1995). 
Generally, these are based on observations at a single epoch. Nevertheless, because the 
ephemerides were well known and because, in any event, the amplitudes of RRab’s in K 
band are only <0.3 mag, Carney et al. (1995) believe that the individual photometry errors 
are only 0.03 mag. 


5. Determination of the Zero Point 

Figure 3 shows the individual estimates for A Ay,i for the 20 stars for which data are 
available. Note that A Ay,i is (1 — a) -1 times the difference between the last two columns 
in Table 1 (see eq. (|2|)). The mean value is (A Ay) = —0.11 with a scatter of 0.22 mag. The 
best estimate of the offset to the Stanek (1996) map is therefore, 

A A v = —0.11 ± 0.05 (RR Lyraes). (5) 

The error in this estimate is determined from the scatter. In principle, there are also 
uncertainties in the zero points of MACHO photometry and of equation (|T). However, 
from the discussions in §3 and §4.1, we conclude that these zero-point errors are negligible 
compared to equation (j5j). 

As discussed in §3, it is also possible to incorporate metallicity information when 
predicting (V — K ) colors. If we do so, however, we find that we recover equation (|5|) 
exactly. That is, metallicity contains no additional information over and above that already 
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contained in the period. This may well be because the errors in the Walker & Terndrup 
(1991) metallicity measurements are larger than the scatter in the period-color relation 
shown in Figure 1. In any event, since metallicity information does not improve (or change) 
the determination, we ignore it. 

The scatter in Figure 3 can be roughly accounted for as follows. There is a scatter 
of 0.06 mag in ( V — K) from the calibrating relation ([||). The Stanek (1996) map has an 
error of 0.10(1 — a) = 0.09 in E(V — K ). The error in K and V photometry are 0.03 mag 
and 0.12 mag, respectively (see §4). Adding these values in quadrature and multiplying 
by (1 — a) -1 yields a predicted scatter of 0.18 mag in A Ay,i compared to 0.22 ± 0.04 mag 
actually observed. 


6. Discussion 

Equation (H) based on RR Lyrae stars (type ab) agrees very well with equation (|3]) 
based K giants. This is heartening because the largest discrepancy among all methods of 
determining E(B — V)q was between RR Lyraes, E(B — V)o = 0.60 ± 0.03 derived by 
Walker & Mack (1986), and cool giants, E(B — V ) 0 = 0.45 ± 0.04 derived by van den 
Bergh (1971). The same discrepancy is not present when similar classes of stars are used to 
determine the extinction via measurement of E(V — K). We conclude that the E(V — K ) 
method has passed an important test for systematic errors. We therefore combine the two 
determinations to obtain, 


A Ay = —0.11 ± 0.04 (combined), (6) 

as our best overall estimate. This corresponds to Ay = 1.36 for the “Blanco A region” 
originally used by Stanek (1996) to calibrate his map. 

Our approach does not permit a direct measurement of E(B — V)q for this region. 
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However, if we follow Terndrup et al. (1995) and adopt A V /E(B — V) 0 = 2.85, we derive 
E(B — V)o ~ 0.48 which is closer to the value derived by van den Bergh (1971) from cool 
giants than it is to the value derived by Walker & Mack (1986) from RR Lyrae stars. 

We are very grateful for the skilled support by the technical staff at Mount Stromlo 
Observatory. Work at LLNL is supported by DOE contract W7405-ENG-48. Work at the 
CfPA is supported NSF AST-8809616 and AST-9120005. Work at MSSSO is supported 
by the Australian Department of Industry, Technology and Regional Development. Work 
at OSU is supported in part by grant AST 94-20746 from the NSF. WJS is supported by 
a PPARC Advanced Fellowship. KG thanks support from DOE OJI, Sloan, and Cottrell 
awards. CWS thanks support from the Sloan, Packard and Seaver Foundations. 
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Table 1: Columns 1 and 2 show the positions of each star in 2000 coordinates. Column 3 
gives the period and column 4 gives the apparent V magnitudes, both from MACHO data. 
Column 5 gives the apparent K magnitudes from Carney et al. (1995). Column 6 gives the 
metalicity from Walker & Terndrup (1991). Column 7 gives the visual extinction from the 
Stanek (1996) map. Columns 8 and 9 gives the actual (V — K) (from columns 4,5 and 7) 
and the (V — K ) predicted. 
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Fig. 1.— Period-color relation for 20 nearby RR Lyrae stars including 17 RRab’s ( circles ) 
and 3 RRc’s {crosses). Data are taken from Jones et al. (1992). Straight line is the best fit 
to the RRab’s: (V — K) 0 = 1.407 + 1.245 logPo- Note that RRc’s do not follow this relation. 
They are therefore excluded from further consideration. 

Fig. 2.— Difference between MACHO and OGLE photometry {Vmacho ~ Vogle) for 50 
RRab’s as a function of MACHO magnitude, Vmacho ■ The OGLE data are from Olech 
(1997) and have been corrected from mean magnitude to magnitude at mean flux. The 
mean offset is ( Vmacho ~ Vogle) = 0.006 ± 0.019. 

Fig. 3.— Individual estimates of the zero-point offset to the Stanek (1996) map, A Ay^ for 
20 RRab’s. See Table 1 and equation (Q). The mean offset is A Ay = —0.11 ± 0.05. 
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